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^ ! Abstract 

C^ , We present the main results of our microcanonical simulation of the 

I ' Wess Zumino Witten action functional. This action, being highly non- 

^~^' trivial and capable of exhibiting many different phase transitions, is 

r~| , chosen to be representative of general complex actions. We verify the 

^ I applicability of microcanonical simulation by successfully obtaining 

two of the many critical points of the Wess Zumino Witten action. The 

microcanonical algorithm has the additional advantage of exhibiting 
a I critical behaviour for a small 8x8 lattice. We also briefly discuss 

the subtleties that, in general, arise in simulating a complex action. 

Our algorithm for complex actions can be extended to the study of 

D-branes in the Wess Zumino Witten action. 

1 Description and Results 

In Minkowski spacetime, the Feynman path integral has a complex integrand 
given by e* *^ , where Sm is the Minkowski action. On Euclidean continuation 
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of spacetime, the integrand becomes real, namely, e~^. The positive definite 
Euclidean action Se is ideally suited for Monte Carlo simulations based on the 
canonical ensemble. There is, however, a class of actions that are generically 
of the form S = Sr + iSj in both Minkowski and Euclidean spacetime. These 
actions are of fundamental importance in physics, and we refer to them as 
complex actions. For example, all topological effects appear in the action in 
the form of complex actions. A chemical potential term renders the action 
complex. Lattice QCD at finite temperatures and and densities has a complex 
action, and so do chiral gauge theories. Complex actions appear in study of 
non-Hermitian hamiltonians, complex Langevin processes, random matrix 
models and so on. 

A complex valued action does not have any probabilistic interpretation, 
and hence, for example the Metropolis algorithm is unsuitable. All attempts 
to numerically simulate complex actions for quantum field theories have met 
with formidable difficulties |l|, H, |^ due to diverging errors arising from in- 
creasingly large volumes needed for obtaining the continuum limit. In partic- 
ular, the conventional Metropolis algorithm fails for complex actions due to 
large errors in simulating the expectation values of cosines and sines of exten- 
sive quantities. These expectation values are necessary due to the presence 
of the imaginary component in the action. 

There have been some studies that simulate complex actions using the 
Langevin equation; however, the results of these simulations cannot be shown 
to converge to the required distribution g, |^ . Other methods ^ also suffer 
from either a restriction to a very small lattice, or are not efficient or well 
defined. 

In our work, we work directly with the Feynman path integral formulation 
of quantum field theory. We apply a generalized form of the microcanoni- 
cal method to the case of complex actions. The microcanonical method 
has been very successful in simulating global quantities of real statistical 
systems [0, H, |^. The microcanonical method yields a direct computation 
of the partition function (up to a proportionality constant) and hence the 
determination of its (complex) zeros. The zeros themselves contain useful 
information, in particular whether the system has critical points. The mi- 
crocaconical simulation technique therefore offers an attractive alternative 
to the usual methods, based on the canonical ensemble, for determining the 
phase diagram of statistical systems. 

In this paper we extend the microcanonical algorithm to complex actions, 
and, compute its partition function. In particular, we study the well known 



and highly nontrivial Wess Zumino Witten (WZW) complex valued action 
functional. Since this model exhibits many distinct phase transitions, we 
are able to provide a solid test of our algorithm. Our simulations predicts 
the critical parameters, and we confirm the effectiveness and accuracy of our 
simulation by comparing the simulated results with the known analytical 
results. 

The Wess Zumino Witten action in two-dimensional Euclidean space has 
the following form 

S = Sx + ikSi, (1) 

where the chiral term is given by 

Sx = ^Jd'xTr{d,gd,g-') (2) 

and the imaginary piece comes from the Wess Zumino term 

24:nJ ^ ^ ^^ dy'^ dy^^ dy^^ ^ ' 

S3 

For the quantum field theory to be well defined, k must be integer-valued. 

For the case of SU{2), the Wess Zumino term can be parametrized in 
terms of the Euler angles [|T^] , namely, 

Si = - f £x(j){x)sin^'ilj{x)sine{x)e'"'dMx)dJ{x). (4) 

n J 

As is well known, the WZW model has infinitely many critical points, 
one for each integer value for k. The WZW model is conformally invariant 
at the critial coupling constants given by[|10 
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Ag = — ; k : integer. (5) 

By extending the Lee- Yang theorem for zero's of the partition function 
to the complex case, we expect that the partition function, in the infinite 
volume limit, will be zero close to the critical coupling constants. We should 
be able to locate a zero (close to the continuum model) if the size of the 
lattice is large enough. 

By calculating the simulated partition function as a function of k, the 
zeroes can be identified as the points of the k axis where the partition function 



Z crosses over from a positive to a negative value. This cross over is a unique 
feature of complex actions, since the partition function of a real action can 
never be negative. Given that only integer fc's are physical, the only physical 
zeros of Z are those that are close to integer values of k. 

The partition function of the Wess Zumino Witten action functional is 
given by 

Z(A, A;) = /"[£>[/] exp(-5A) exp(zfc5/), (6) 

which in turn can be cast into the form 

/ + 00 
dEp{\ E) exp{ikE) (7) 

-oo 

with the density of states given by 

p(A, E) = I [DU]6iE - Sj) exp{-S,) (8) 

The density of states is estimated by a generalization of the microcanon- 
ical algorithm. First, the simulation space is broken up into energy bins 
of size AE. Bin E = is taken to be the reference bin. The system is 
thermalized by randomly generating configurations and accepting only those 
configurations whose values of 5*/ fall within the first two bins. The density 
of states p(A, AE) is estimated by calculating the occupation number ratio 
between the two bins. The updating of configurations and crediting of counts 
to the bins are based on a combination of the microcanonical and canonical 
ensembles. This is done as follows. 

1. A new configuration is generated, and its acceptance is weighted by the 
usual Metropolis factor exp(— S'a). Note that since the chiral contribu- 
tion, being real, plays the role of a weighting action, the value of A has 
to be at the critical value for a particular level; otherwise, if A 7^ Ac, 
the system would never go critical. 

2. The updated configuration is accepted only if it's Si value falls within 
the two bins. If this is so, the counter for the bin where Sj lands is 
credited. 

3. If the updated configuration does not satisfy the above criterion, a 
count will still be credited to the bin of the old configuration. 




Figure 1: Density of states for two values Ac 



By continuing this process across all the energy bins allows us to get an 
estimation of the density of states. A crucial point in going from one set of 
bins to another is to have, in between, at least one overlapping bin for the 
purpose of normalizing the occupation number of the various bins. 

All of the results quoted are for simulations for a 8 x 8 lattice, as this 
size is sufficient for accurately estimating the value of critical k. In our work, 
we tested our simulations for two levels, namely k = 1 and k = 2. For the 
critical coupling A^ = An,AE was set to be 0.2 while a finer value of 0.1 
was chosen for A^ = 2n. A total number of about 1 million configurations 
were generated for every set of energy bins. The energy bin was terminated 
when the resulting density of state is less than 3 decimal places (which is the 
achieved level of accuracy in our simulations). The results for the simulated 
density of states is shown in Figure |I[ 

Since p{X,—E) = p{X,E), it is only necessary to perform the simula- 
tion for E > 0. This refiection symmetry ensures that the WZW partition 
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Figure 2: The integrand of Z for Ac = 47r 



function is real, since it simplifies to 



Z(A, k) = 2 dEp{X, E)cos{kE) 





(9) 



The simulated value for integrand of Z{Xc, k) is given in Figures |] and |^. 
The plot of Z{Xc, k) against k is physical relevant only for integer fc's. We 

expect that Z = at the value of k that corresponds to A^ = -;-; for other 

k 
values of integer k, Z is physical and contains information on how the theory 

behaves away from criticality. 

The behaviour of Z near its zero is given in Figures § and ^. By scanning 

for zeros of the simulated partition function Z that are close to integers, we 

find that: 

A^simuiated = 1-06 ± 0.07 (for A^ = 47r) 

and 

A:simuiated = 1-95 ± 0.06 (for A^ = 2^). 
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Figure 3: The integrand of Z for Ag = lix 



The errors were deduced from the standard deviation computed from two 
separate sets of results for each couphng. 

The exact values are k = \ and k = 1 respectively, and we see that 
we have obtained a result of fairly high accuracy with a surprisingly small 
lattice. 



2 Discussion 



Our simulations of the Wess Zumino Witten action demonstrates that the 
microcanonical method is in principle capable of simulating complex actions. 
Unlike the case of real actions, however, there are several subtleties involved. 
Chief among them are the following. 

1. A high degree of accuracy is required in order to achieve stable and 
consistent results. The reason being that, unlike real actions, terms 
in the complex action partition function cancel each other (see the 
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Figure 4: Scan of Z near a zero for Ac = 47r. The zero occurs at /c ^ 1.08. 



graphs of Figures |^ and |^), leaving behind a small value in the region 
where the phase transitions take place. 

The fact that the critical point can be made to move from /c = 1 to 
A; = 2 by changing the value of Ac in the weighting function is strong 
evidence that our simulation is giving the correct result. There are 
infinitely many critical points of the WZW model, one for each integer 
/c, and in principal, all these critical points can be obtained by our 
simulation. 



3. There is an additional necessity of having to provide a good resolution of 
the cos{kE) term. Typically, this means that there should be at least a 
few points within each trigonometrical quadrant. Consequently, since 

AE ~ — , a smaller energy spacing is required for larger couplings. 

k 
Computational time therefore is expected to increase as coupling k 

increases. 




Figure 5: Scan of Z near a zero for Ac = 27r. The zero occurs at /c ~ 1.89. 



4. Since the cosine of an argument cannot be taken as a rehable weight 
in any Monte Carlo simulation, the contribution to the weight must 
always come from the real action. This constraint demands that the 
coupling A in the real action has to be fixed at its critical value Ac. 
In other words, it is not possible to fix the coupling in the imaginary 
component, namely A;, and deduce the critical coupling Ac of the real 
action. 



Our method is ideally suited for studying the critical properties of 
extended systems. In particular, D-branes in the WZW model are non- 
local background configurations that respect the conformal invariance 



of the theory [|r^. Our method, for example, can numerically study 
under what circumstances a D-brane background in the WZW model 
results in a conformally invariant system. 

In conclusion, we believe that the ability of the microcanonical technique 
to get accurate predictions from the Wess Zumino Witten action is a major 



step forward in the simulation of complex actions, and opens the way for 
numerically studying other complex actions, as well as extended objects such 
as D-branes. 
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